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In the standard model (SM), semileptonic decays of B 
mesons are well-understood processes mediated by a W 
boson [H-Q. Decays involving the higher mass r lepton 
are sensitive to additional amplitudes, such as those in- 
volving an intermediate charged Higgs boson [iHII, and 
offer an excellent opportunity to search for this and other 
non-SM contributions. 

Our understanding of exclusive semileptonic decays 
has greatly improved over the past two decades, thanks 
to the development of heavy-quark effective theory and 
precise measurements of B — » Dt*H~T>i || at the B fac- 
tories Hlj |. SM expectations for the relative rates 
ft(£>«) = B(B D^t~v t )/B{B -> D^H-VA have 
less than 6% uncertainty Calculations 041] based 
on two-Higgs-doublet models (2HDM) predict a substan- 
tial impact on the ratio 7^(1?) and a smaller effect on 
K(D*). The ratios K(D) and TZ(D*) are independent of 
the Cabibbo-Kobayashi-Maskawa element \V c b\ and also, 



to a large extent, of the parameterization of the hadronic 
matrix elements. 

The decay B — > D*t~V t was first observed in 2007 
by the Belle Collaboration [l2| . Since then, both BABAR 
and Belle have published improved measurements and 
have found evidence for B — > Dt~T> t decays [13415) . Al- 
though the measured values for 1Z(D) and 1Z(D*) have 
consistently exceeded the SM expectations, the signifi- 
cance of the excess has remained low due to the large 
statistical uncertainties. 

This analysis is an update of an earlier BABAR mea- 
surement [l3| . It is based on the full BABAR data sample 
and includes improvements to the event reconstruction 
that increase the signal efficiency by more than a factor 
of 3. 

We analyze data recorded with the BABAR detector [nl 
at a center-of-mass (cm.) energy of 10.58 GeV, cor- 
responding to the mass of the T(AS) resonance, which 
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decays almost exclusively to BB pairs. The data sam- 
ple comprises an integrated luminosity of 4261b" 1 , and 
contains 471 x 10 6 BB pairs. An additional sample of 
40fb _1 , taken at a cm. energy 40 MeV below the T(45) 
resonance (off-peak data), is used to study continuum 
background from e + e~ — > //(7) pair production with 
/ = u,d,s,c,T. 

We choose to reconstruct only the purely leptonic de- 
cays of the t lepton, r~ — > e~V e v T and t~ — > /x _ z7„i/ r , 
so that B -> D^t~v t (signal) and B -> D^H~Vi 
(normalization) events are identified by the same par- 
ticles in the final state. This leads to the cancellation 
of various sources of uncertainty in the ratios TZ(D^). 
Events corresponding to T(45) — > BB decays are se- 
lected by reconstructing the hadronic decay of one of 
the B mesons (B tag ), a meson and a lepton (e 

or ji). Signal and normalization yields are extracted 
from a fit to the spectra of two variables: the invari- 
ant mass of the undetected particles m miss = p miss = 
(Pe+e- ~ Ptag — PdM ~ Pi) 2 (where p denotes the four- 
momenta of the colliding beams, the B tag , the D^*\ and 
the charged lepton) and the lepton three-momentum in 
the B rest frame \p* t \ . The m miss distribution of decays 
with a single missing neutrino peaks at zero, whereas 
signal events, which have three missing neutrinos, have a 
broad ro miss distribution that extends to about 9GcV 2 . 
The observed lepton in signal events is a secondary par- 
ticle from the r decay, so its |p| | spectrum is softer than 
for normalization events. 

The Stag reconstruction has been greatly improved 
with respect to previous analyses fl7| . We now recon- 
struct B tag candidates in 1680 final states. We look for 
decays of the type Stag — > SX^, where S refers to a seed 
meson (D°, D*°, D+, D*+ , D+, D*+, or J/ip) recon- 
structed in 56 different decay modes, and is a charged 
state decaying to up to five hadrons (tt ± , K ± , 7T°, and 
Kg). Two kinematic variables are used to select B tag can- 
didates: vies '- 



F 2 

beam 



P?a 



and AE = E, 



tag ^bcam- 



Here p tag and E^ ag refer to the cm. momentum and en- 
ergy of the -Btag, and £?beam is the cm. energy of a single 
beam particle. For correctly reconstructed B decays, the 
wies distribution is centered at the B-meson mass with a 
resolution of 2.5 MeV, while AE is centered at zero with 
a resolution of 18 MeV. We require tues > 5.27 GeV and 
\AE\ < 0.072 GeV. 

We combine each B tag candidate with a D^*' meson 
candidate and a charged lepton I. Events with addi- 
tional charged particles are rejected. The laboratory 
momentum of the electron or muon is required to 
exceed 300 MeV or 200 MeV, respectively. D decays 
are reconstructed in the following decay modes: D° — > 
K-ir + J<-K+ 1 K-ir+ir° 1 K-Tr + ir-ir + , K° s tt+tt-- D+ -> 

R-TT+TT+, K-TT+TT+IT , A"°7T+ , TT+ Tr~ , K° S 1T + TT° , 

K°K+, with K° s -> tt+tt". For D* candidates, the 
decays D*+ -> D°ir+, D+ir°, and D*° -> D°n ,D°j are 



used. 

In events with more than one reconstructed BB pair, 
we select the candidate with the lowest value of E ex $ Ia , 
defined as the sum of the energies of all photon candidates 
not associated with the reconstructed BB pair. We fur- 
ther reject combinatorial background and normalization 
events by requiring q 2 = (pb — P_d<»)) 2 > 4 GeV 2 and 
|p miss | > 200 MeV, where |p m j ss | is the missing momen- 
tum in the cm. frame. 

We divide the B tag D^*>£ candidates that satisfy the 
previous requirements into eight subsamples: four D^*'£ 
samples, one for each of the types of charm meson recon- 
structed (D°£, D*°£, D+£, and D*+£), and four D^tt°£ 
control samples with the same selection plus an addi- 
tional 7T°. These control samples constrain the poorly un- 
derstood B — > D**(1/t)v background (where D** refers 
to charm resonances heavier than the D* ground state 
mesons), which enters the D^£ sample predominantly 
when the ir° from D** — > D^^ir decays is not de- 
tected. The D^t: £ samples have a very large contin- 
uum background, so we restrict this sample to events 
with | cos A6*thrust| < 0.8, where A6* t hrust is the angle be- 
tween the thrust axes of the B tag and of the rest of the 
event. 

We improve the separation between well-reconstructed 
events (signal and normalization) and the various back- 
grounds by using boosted decision tree (BDT) selectors 
[181 ] . For each of the four D^£ samples, we train a 
BDT to select signal and normalization events and re- 
ject D**£U background and charge cross-feed, defined as 
D^*\£/t)v decays reconstructed with the wrong charge. 
Each BDT selector relies on the simulated distributions 
of the following variables: (1) B ox t ra ; (2) AE; (3) the re- 
constructed mass of the signal D meson; (4) the mass dif- 
ference for the reconstructed signal D*: Am = m(Dir) — 
m(D); (5) the reconstructed mass of the seed meson of 
the -Btag; (6) the mass difference for a D* originating 
from the B tag , Am tag = m(-D tag 7r) - m(D tag ); (7) the 
charged particle multiplicity of the B tag candidate; and 
(8) cos Aflthrust- For the D^tt°£ samples, we use simi- 
lar BDT selectors that are trained to reject continuum, 
DW(£/t)u, and other BB events. After the BDT re- 
quirements are applied, the fraction of events attributed 
to signal in the m miss > 1.5 GeV 2 region, which excludes 
most of the normalization decays, increases from 2% to 
39%. The remaining background is composed of nor- 
malization events (10%), continuum (19%), D**(£jr)v 
events (13%), and other BB events (19%), primarily from 
B -> D^D ( * )+ decays with D+ t+v t . 

As described below, the fit procedure relies on the 
Monte Carlo (MC) simulation [lj|-[2l[ of the two- 
dimensional m^^-lp} | spectra of the different signal and 
background contributions. For semileptonic decays, we 
parameterize the hadronic matrix elements of the sig- 
nal and normalization decays by using heavy-quark ef- 
fective theory-based form factors (FFs) 22[. For low- 
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mass leptons, there is effectively one FF for B — > Di~vi 1 
whereas there are three FFs for B — > D*i~vg decays, all 
of which have been measured with good precision [23|. 
For heavy leptons, each of these decays depends on an 
additional FF which can be calculated by using heavy- 
quark symmetry relations or lattice QCD. We use the 
calculations in Ref. Q for B — > Dt~V t and in Ref. Q 
for B — > D*t~v t . For the D**(£/t)v background, we 
consider in the nominal fit only the four L = 1 states 
that have been measured [HI ■ We simulate these decays 
by^ using the Leibovich-Ligeti-Stcwart-Wise calculation 

We validate and, when appropriate, correct the sim- 
ulations by using three data control samples selected 
by one of the following criteria: Ecxtra > 0.5 GeV p6| . 
q 2 < 4 GeV 2 , or 5.20 < m E s < 5.26 GeV. We use off- 
peak data to correct the efficiency and the spectrum 
of simulated continuum events. After this correction, the 
m^jgg and \p}\ distributions of the background and nor- 
malization events agree very well with the simulation. 
However, we find that small differences in the £^ G xtra spec- 
trum and other BDT input distributions result in a 5%~ 
10% efficiency difference between selected data and MC 
samples. We correct the continuum and BB backgrounds 
by using the 5.20 < wes < 5.26 GeV control sample. The 
same correction, with larger uncertainties, is applied to 
D**{£/t)v events, since their simulated Eextra spectrum 
is very similar. 

We extract the signal and normalization yields from 
an extended, unbinned maximum-likelihood fit to two- 
dimensional TimissHpH distributions. The fit is per- 
formed simultaneously to the four D^*'£ samples and the 
four D^TT^tt samples. The distribution of each D^H 
sample is described as the sum of eight contributions: 
Dtv, D*tv, Dlv, D*£v, D**(1/t)v, charge cross-feed, 
other BB, and continuum. The yields for the last three 
backgrounds are fixed in the fit to the expected values. 
A large fraction of B — > D*iv decays (for B = B° or 
B + ) is reconstructed in the D£ samples (feed-down). We 
leave those two contributions free in the fit and use the 
fitted yields to estimate the feed-down rate of B —> D*tv 
decays. Since B — > D{£/t)v decays contributing to the 
D*£ samples are rare, their rate is fixed to the expected 
value. 

The four D^tt°£ samples are described by six contri- 
butions: The D^tv and D^tv yields are combined but 
otherwise the same contributions that describe the D^*'£ 
samples are employed. The four D**{1/t)v yields in the 
control samples are free in the fit, but their ratios to the 
corresponding D**(£/t)v yields in the D^*H samples are 
constrained to the expected values. 

The fit relies on 8x4 + 6x4 = 56 probability density 
functions (PDFs), which are determined from MC sam- 
ples of continuum and BB events equivalent to 2 and 9 
times the size of the data sample, respectively. The two- 
dimensional to 2 ^ s -|p|| distributions are described by us- 




FIG. 1. (Color online) Comparison of the data and the fit 
projections for the four D^t samples. The insets show the 
\P*t \ projections for m^iss > lGeV 2 , which excludes most of 
the normalization modes. In the background component, the 
region above the dashed line corresponds to charge cross-feed, 
and the region below corresponds to continuum and BB. 



ing smooth nonparametric kernel estimators |27( ■ The fit 
is iterated to update some of the parameters that depend 
on the normalization yields, most importantly the rate of 
signal feed-down. This process converges after the first 
iteration. We performed MC studies to verify that nei- 
ther the fit procedure nor the PDFs produced significant 
biases in the results. 
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Figure Q] shows the and projections of the fit 

to the four D^£ samples. The fit describes the data well 
and the observed differences are consistent with the sta- 
tistical and systematic uncertainties on the signal PDFs 
and background distributions. 

We extract the branching fraction ratios as 1Z(D^) = 
(N s i s /N noIm )/(e sig /£ nolm ), where AT si g and N norm refer 
to the number of signal and normalization events, re- 
spectively, and e S ig/£norm is the ratio of their efficiencies 
derived from simulations. Table U shows the results of the 
fits for the four individual samples as well as an additional 
fit in which we impose the isospin relations 1Z(D°) = 
K{D+) = 11(D) and 1Z(D*°) = 1Z(D*+) = K(D*). The 
statistical correlations arc —0.59 for 1Z(D°) and 1Z(D*°), 
-0.23 for K{D+) and K(D*+), and -0.45 for 11(D) and 
1Z(D*). We have verified that the values of IZ(D^) from 
fits to samples corresponding to different run periods are 
consistent. We repeated the analysis varying the selec- 
tion criteria over a wide range corresponding to changes 
in the signal-to-background ratios between 0.3 and 1.3, 
and also arrive at consistent values of IZ(D^). 

The largest systematic uncertainty affecting the fit re- 
sults is due to the poorly understood B — > D**(1/t)v 
background. The PDFs that describe this contribution 
are impacted by the uncertainty on the branching frac- 
tions of the four B — >• D**lv decays, the relative 7r°/7r ± 
efficiency, and the branching fraction ratio of B — > D**tv 
to B — > D**£v decays. These effects contribute to an un- 
certainty of 2.1% on 11(D) and 1.8% on 1Z(D*). We also 
repeated the fit including an additional floating compo- 
nent with the distributions of B — > D^rjiv, nonresonant 
B ->■ D^Tr(-ir)£v, and B -)■ D**(-^ Tm)tv decays. 
The B — > D**(1/t)v background is tightly constrained 
by the D^tt°£ samples, and, as a result, all these fits 
yield similar values for IZ(D^). We assign the observed 
variation as a systematic uncertainty, 2.1% for 1Z(D) and 
2.6% for 1Z(D*). 

We also account for the impact of the uncertainties 
described above on the relative efficiency of the B —> 
D**(£/t)v contributions to the signal and D^tt°£ sam- 
ples. In addition, the BDT selection introduces an un- 
certainty that we estimate as 100% of the efficiency cor- 
rection that we determined from control samples. These 
effects result in uncertainties of 5.0% and 2.0% on H(D) 
and 1Z(D*), respectively. 

The largest remaining uncertainties are due to the con- 
tinuum and BB backgrounds [4.9% on 1Z(D) and 2.7% on 
1Z(D*)], and the PDFs for the signal and normalization 
decays (4.3% and 2.1%). The uncertainties in the effi- 
ciency ratios £ S ig/£norm are 2.6% and 1.6%; they do not 
affect the significance of the signal and are dominated by 
the limited size of the MC samples. Uncertainties due 
to the FFs, particle identification, final-state radiation, 
soft-pion reconstruction, and others related to the detec- 
tor performance largely cancel in the ratio, contributing 
only about 1%. The individual systematic uncertainties 
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FIG. 2. (Color online) Comparison of the results of this anal- 
ysis (light gray, blue) with predictions that include a charged 
Higgs boson of type II 2HDM (dark gray, red). The SM cor- 
responds to tan/3 /m H + = 0. 



are added in quadrature to define the total systematic 
uncertainty, reported in Table [I] 

There is a positive correlation between some of the 
systematic uncertainties on 1Z(D) and 1Z(D*), and, as a 
result the correlation of the total uncertainties is reduced 
to -0.48 for 1Z(D°) and 1Z(D*°), to -0.15 for 1Z(D+) and 
K(D*+), and to -0.27 for 11(D) and 1Z(D*). 

The statistical significance of the signal is determined 
as Est&t = -\/2A(hbC), where A(ln£) is the change in the 
log-likelihood between the nominal fit and the no-signal 
hypothesis. The statistical and dominant systematic un- 
certainties are Gaussian. We estimate the overall signifi- 

•^stat X CTstat/ \f®s 



cance as St, 



^tot = ^stat A "stat/ Y CTstat 1" "syst! wuere G> sta t 

is the statistical uncertainty and crg yst is the total system- 
atic uncertainty affecting the fit. The significance of the 
B — > Dt~T> t signal is 6. 8c, the first such measurement 
exceeding 5er. 

To compare the measured 1Z(D^*>) with the SM pre- 
dictions we have updated the calculations in Refs. [H |3l[ 
taking into account recent FF measurements. Averaged 
over electrons and muons, we find 1Z(D)sm = 0.297 ± 
0.017 and K(D*) SM = 0.252 ±0.003. At this level of pre- 
cision, additional uncertainties could contribute 0], but 
the experimental uncertainties are expected to dominate. 

Our measurements exceed the SM predictions for 
11(D) and K(D*) by 2.0cr and 2.7cr, respectively. The 
combination of these results, including their —0.27 cor- 
relation, yields \ 2 = 14.6 for 2 degrees of freedom, corre- 
sponding to a p value of 6.9 x 10~ 4 . Thus, the possibility 
of both the measured H(D) and 1Z(D*) agreeing with the 
SM predictions is excluded at the 3.4er level |3_2j. 

Figure [2] shows the effect that a charged Higgs bo- 
son of the type II 2HDM @, would have on 11(D) 
and 1Z(D*) in terms of the ratio of the vacuum expecta- 



where <j s 
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TABLE I. Results of the isospin-unconstrained (top four rows) and isospin-constrained fits (last two rows). The columns 
show the signal and normalization yields, the ratio of their efficiencies, TZ(D^), branching fractions, and E s t&t and E to t, the 
statistical and total significances, respectively. Where two uncertainties are given, the first is statistical and the second is 
systematic. The branching fractions B(B — > D^*\~V T ) are calculated as ft(L> w ) x B(B -> D^Tue), using the average 
B —J- D^t~Vt branching fractions measured by BABAR [2^ - f30( ] . The stated branching fractions for the isospin-constrained fit 
refer to B~ decays. 



Decay 






norm 


£sig/ £norm 


ft(D w ) 


B(B -> D M ru)(%) S atat 


■£tot 


B-^ 




314 ± 60 


1995 ± 55 


0.367 ± 0.011 


0.429 ± 0.082 ± 0.052 


0.99 ± 0.19 ± 0.13 


5.5 


4.7 


B~ — ► 


D*°t-V t 


639 ± 62 


8766 ± 104 


0.227 ± 0.004 


0.322 ± 0.032 ± 0.022 


1.71 ± 0.17 ± 0.13 


11.3 


9.4 


B° -s- 




177 ± 31 


986 ± 35 


0.384 ± 0.014 


0.469 ± 0.084 ± 0.053 


1.01 ± 0.18 ± 0.12 


6.1 


5.2 


B° -»■ 


D* + T~V T 


245 ± 27 


3186 ± 61 


0.217 ± 0.005 


0.355 ± 0.039 ± 0.021 


1.74 ± 0.19 ± 0.12 


11.6 


10.4 


B ->■ 


Dt~V t 


489 ± 63 


2981 ± 65 


0.372 ± 0.010 


0.440 ± 0.058 ± 0.042 


1.02 ± 0.13 ± 0.11 


8.4 


6.8 


B -> 


D*T~V T 


888 ± 63 


11953 ± 122 


0.224 ± 0.004 


0.332 ± 0.024 ± 0.018 


1.76 ± 0.13 ± 0.12 


16.4 


13.2 



tion values tan/3 = V^/vi, and the mass of the charged 
Higgs mjj+ . We estimate the effect of the 2HDM on our 
measurements by reweighting the simulated events at the 
matrix element level for 20 values of tan/3/m^+ over the 
[0.05, 1] GeV -1 range. We then repeat the fit with up- 
dated PDF shapes and £ S i g /£norm values. The increase 
in the uncertainty on the PDFs and the efficiency ra- 
tio is estimated for each value of tan/3/m#+. The other 
sources of systematic uncertainty are kept constant in 
relative terms. 

The measured values of TZ(D) and 1Z(D*) match 
the predictions of this particular Higgs model for 
tanj3/m H + = 0.44±0.02 GeV -1 and tan/3 /m H+ = 0.75± 
0.04 GeV -1 , respectively. However, the combination of 
11(D) and TZ(D*) excludes the type II 2HDM charged 
Higgs boson with a 99.8% confidence level for any value 
of tan/3/m^+. This calculation is valid only for values 
of rriH+ greater than about 10 GeV 0, 0]- The region for 
m H + < 10 GeV has already been excluded by B — > X s j 
measurements [IH, and, therefore, the type II 2HDM is 
excluded in the full tan/3-m ff + parameter space. 

In summary, we have measured the B — > Dt~~U t and 
B — } D*t~v t decays relative to the decays to light lep- 
tons B -> D(*h~V e . We find 

K{D) = 0.440 ± 0.058 ± 0.042 
K{D*) = 0.332 ± 0.024 ± 0.018 . 

These results supersede the previous BABAR results and 
have significantly reduced uncertainties. The measured 
values are compatible with those measured by the Belle 
Collaboration |!l0[l5|. 

The results presented here disagree with the SM at 
the 3.4cr level, which, together with the measurements 
by the Belle Collaboration, could be an indication of 
new physics processes affecting B —> D^*'t~V t decays. 
However, our results are not compatible with the widely 
discussed type II 2HDM for any value of tan/3 and m H + . 
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